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An experimental analysis of the vortex-induced vibrations of a hanging string with variable tension 
along its length is presented in this paper. It is shown that standing waves develop along the hanging 
string. The evolution of the Strouhal number St with the Reynolds number Re first follows a trend 
similar to what is observed for a circular cylinder in a flow for relatively low Reynolds numbers 
(32 < Re < 700). Second, the extracted mode shapes are self-similar : a rescaling of the spanwise 

11^ \ coordinate by a self-similarity coefficient allows all of them to collapse on a unique function. The 

self-similar behaviour of the spatial distribution of the vibrations along the hanging string is then 
explained theoretically by performing a linear stability analysis of an adapted wake-oscillator model. 

^^ ^ This linear stability analysis finally provides an accurate description of the mode shapes and of the 

'^H , evolution of the self-similarity coefficient with the flow speed. 

Ch ■ I. INTRODUCTION 

I ' Vortex-induced vibrations (VIV) of slender structures are of major importance in many engineering fields. In 

^ , offshore applications for instance, a good knowledge of the dynamic loads exerted on risers and mooring cables is 

t^_H ' required to prevent fatigue damage [l|, y] • VIV are self-sustained oscillations of an immersed bluff body originating 

c/3 in a strong coupling between the structure's dynamics and its induced fluctuating wake. This coupling may lead to a 

O ' synchronization of the solid vibrations and vortex shedding frequencies, a phenomenon known as lock-in, resulting in 

t^ _ large amplitude oscillations [3, |j| . 

J>~^. The case of long and flexible structures in VIV, most relevant to offshore applications, has been extensively studied 

i-^ ] using numerical computations @, Q, experiments [3, Q and reduced-order models @, UOl- When a long flexible 

I '~^ , structure is placed perpendicular to the oncoming flow, higher vibration modes can be excited as the flow vel ocity 

is increased, through the successive synchronizations of its structural modes with the local shedding frequency [llj |. 

T-H , Spatial variations of the vortex shedding frequency, for example in sheared flows, significantly modifies the dynamics 

K*" ■ from the canonical problem of tensioned cables in uniform flows |12h14J . Similarly, spatial variations of the cable's 

^^ ] tension (for example due to gravity and/or free end conditions) is expected to impact the structure's dynamics, in 

particular by localizing vibrations in regions with lower tension. This second configuration has so far received little 

attention jlSl . Ilq . The present paper thus discusses an experimental investigation of the VIV of a hanging string, 

. where gravity-induced tension does vary from bottom to the top. 

^i Section [IT] presents the experimental setup and typical raw measurements. In Section IIIIl the evolution of the 

iTjr , dimensionless frequency of oscillations with Reynolds number is analyzed and Section IIVI identifies a self-similar 

I ■ behaviour in the spatial distribution of the amplitude of vibrations. This self-similarity is confirmed analytically in 

:• [ Section |V] using a wake-oscillator model, before concluding remarks are given in Section IVll 

X: 

'^lJ . II. EXPERIMENTS 

.9.' 

A. Experimental setup 

Experiments were performed in the Fluid-Structure Interactions laboratory of the University of Massachussets in 
Amherst. Fixed at its top extremity, a chain was hung in the water channel and was submitted to a flow of constant 
speed U, figure [1] (a). The chain consisted of several spheres of equal diameter connected through short metalhc 
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FIG. 1: (a) Experimental setup of the hanging string subject to a transverse flow. (6) Geometry of the chains that were used 
to model a string of equivalent diameter Deq. (c) Visualization of Chain 3 during one typical experiment. 

L (m) D (mm) d (mm) h (mm) nis (g/m) 

Chain 1 0.48 2.0 3.0 0.4 8 

Chain 2 0.48 3.1 4.5 0.6 21 

Chain 3 0.48 4.4 6.3 0.8 29 

TABLE I: Dimensional properties of the three different chains used for experiments. 



rods, figure [T] (6). Three different chains were used, and their respective characteristics are shown in table H] where 
L and D are the chain length and diameter, d is the distance between two successive spheres, h the diameter of the 
connecting metallic rods and rris the chain mass per unit length. The cross-flow displacement Y(Z,T) was recorded 
using a Phantom Miro MHO camera through a window located at the downstream end of the water channel, figure 
[T] (a). The camera was equipped with a lens of focal length 50 mm and aperture f/1.4. The position of the chain 
was obtained by tracking the contrast between the bright spheres and dark background, figure [T](c), and filtering the 
resulting displacement signal around its dominant frequency. Some wake visualizations were also performed using dye 
released slightly upstream from the lower tip of the string. 

The chains used in the experiments were heavy enough not to float while keeping a very small bending stiffness 
for small curvatures. Yet, significant non-linear rigidity may appear for higher curvatures. When the wavelength is 
large in comparison with the spheres spacing d, the chain can be used to model an ideal string of equivalent diameter 
Deq = ttD/4:. Here, D^q is defined such that the equivalent string has the same cross-sectional area as the spheres, 
thereby neglecting the effect of the connecting rods on the fluid/solid coupling. The Reynolds number Re is based on 
this equivalent diameter and reads as 



Re 



UD, 



(1) 



where z^ = 10^^ m^.s^'^ is the kinematic viscosity of water at room temperature . With this definition, the Reynolds 
number varies in these experiments between 32 and 700. 

For all experiments, the flow velocity was sufficiently low for the string's static in-line defiection to be negligible, 
so that the flow may be considered perpendicular to the string. The in-line vibrations were too small to be detected 
in our setup and are therefore not investigated further, although they are known to exist in such configurations [17|. 
The present study consequently focuses on the string's cross-flow oscillations. Finally, it ought to be mentioned here 
that it was not possible to capture the motion of the top 6 cm (of the total 48cm) of the string because of the size 
of the window through which videos were recorded. This does not affect the conclusions and explains why the top 
displacement of the string is not shown later in figure |3] (a) and (b). 







FIG. 2: Wake visualization during one cycle of oscillations of chain 2 for Re = 106. The white dashed line on each picture shows 
the location of the chain at rest, (a) and (c) correspond to the minimum and maximum displacement of the chain respectively. 

B. Wake visualizations 

Wake visualizations were performed in order to confirm that the measured displacements actually are vortex-induced 
vibrations. Four snapshots showing the tip displacement of the hanging string as well as its near wake are shown 
in figure [21 The vortex shedding behind the string is synchronized with the periodic oscillations of the solid: one 
vortex is shed every time the free end of the hanging string reaches its minimum or maximum displacement, as shown 
in figure [21 The frequency of vibrations of the body and the frequency of the vortex shedding are locked, and the 
vibrations result from this strong coupling between the string's dynamics and its induced wake. 

C. String's motion 

The string's transverse displacement was obtained from the videos and is plotted in figures [3 (a) and (b) for two 
typical cases, showing standing waves along the hanging string. The maximum amplitudes are significant - of the 
order of half of the equivalent diameter. The vibrations with the highest amplitudes are concentrated in the lower 
part of the string, where the tension is lower. 

The standing wave can be decomposed into its temporal and spatial components, as shown in figure [3l (c)-(/). As 
the flow velocity increases, a higher structural mode is excited, which corresponds to a higher frequency of oscillations, 
since VIV are due to the lock-in between vortex shedding and one of the string's eigenfrequencies. In the sample cases 
shown in figure [3l the second and the fourth structural modes are excited, respectively, as the flow velocity increases 
from U = 0.033 m.s^^ to U = 0.064 m.s^^. The excitation of higher modes at higher fiow velocities is expected in 
VIV of flexible cylinders [HI, [H. 

III. STROUHAL NUMBER 

The natural frequency of vortex shedding, Fs, is defined here as the frequency at which vortices would be shed if 
the string were fixed and perfectly rigid. The Strouhal number is thus defined as 



F D 

bt-—^, 



(2) 



For flexible structures, VIV orig inate in the lock-in between the frequency of vortex shedding and one of the eigen- 
frequencies of the structure [11] . The measured frequency F of the string's displacement is therefore close, but not 
equal to, the natural frequency of vortex shedding, F^. The Strouhal number, as defined in equation ([2]), may thus 
be estimated here by FDeq/U, and is shown in figure [Has a function of the Reynolds number. 

The Strouhal number St is close to 0.2, which is the value commonly used for cylinders over a wide range of Re 
[3, [23|. Moreover, St increases with Re for Re < 200 before it reaches a constant level slightly higher than 0.2 for 
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FIG. 3; Vortex-induced vibrations of Chain 3 for two typical experiments, U = 0.033 m.s~^ {Re — 115) (left) and U — 
0.064 m.s"'^ {Re = 221) (right). {a),{h) Dimensionless displacement y = Y/Deq (level step : 0.015), {c),{d) Rescaled tip 
displacement, the amplitudes of the temporal components were arbitrarily set to 1, (e),(/) spatial component of the solid's 
displacement, derived as the time-averaged shape of the string over all the time-steps for which the tip displacement is larger 
than 3/4 of the tip r.m.s. amplitude. 

higher values of Re. This trend is very similar to the St vs Re plot that can be found for instance in Ref. ^ or in 
Ref. [2l|. It was also shown in Refs. [22|, 123] that a fit in a scries of l/\/i?e is appropriate to interpolate the variations 
of St in this range of Re. Such a fit is derived for the presented data and is plotted in figure |4l 
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FIG. 4: Strouhal number as a function of the Reynolds number: measurements (symbols), best fit in a series in l/VRe 
(continuous line) and fit introduced in Ref. [23] (dashed line) 



The coefficients proposed in Ref. [2J] also result in a reasonably good fit for the presented measurements, figure S) 
This shows the relevance of this estimation of the Strouhal number as well as the consistency of its variations with Re. 
The chosen definition of D^q may impact the experimental curve of figure |4j therefore explaining part of its deviation 
from the fits derived in Ref. |23l . 

Equation (jH]) is consequently used in the remainder of the paper to estimate the Strouhal number of the discussed 
experimental results. 

IV. MODE SHAPES 
A. Characteristic dimensions of the problem 

The equivalent diameter Deq is a natural characteristic length for the cross-flow displacement Y{Z,T), and the 
natural frequency of vortex shedding ujg — 2-KStU / D^q defines a characteristic frequency for the string's VIV. 

For the spanwise coordinate Z , a characteristic length Zc is derived from a typical phase speed, as follows. The 
tension is induced by gravity g, and a characteristic value for the tension consequently is = rriggZc. The solid is 
completely immersed so its effective mass per unit length reads rrit = rris + pnD^ /A, where an added mass coefficient 

of 1 is assumed, p — 1000 kg.m^ being the fluid density. By defining a phase speed c ~ ^Jm^gZ^Jmt = ujgZc, a 
characteristic length is derived as 



Zc 
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(4) 



The dimensionless set of parameters consequently is 



Y Z 



(5) 



With this scaling, the problem is characterized by one dimensionless coefficient, the string's reduced length t — L/Zc, 
which depends on the flow speed and thus changes from one experiment to another. 



B. Self-similar mode shapes 

For moderate values of Re, Re < Rec with Rec ~ 110 for Chain 1, Rec — 160 for Chain 2 and Rec ~ 250 for Chain 
3, vortex-induced standing waves are observed. Mode shapes are extracted for every experiment and are normahzed 
using the L^ norm. Figures [3] (e) and (/) show some similarity in the mode shapes of different experiments: the mode 
shape of the lower flow velocity case, figure |3] (e), and the bottom part of the high velocity one {Z < 0.10 m), figure [3] 
(/), actually look alike. We therefore search here for a self-similar behaviour of the spatial distribution of the string's 
vibrations. Considering the mode shape of the experiment with the highest value of ^ as a reference fref (z), the 
spanwise coordinate z is rescaled for each experiment as z/k, so as to minimize the L^ norm of ip {z/ k) — fref {z)'- 

E'^ (k) = / (p (-] - (Pref (z) dz. (6) 

All the rescaled mode shapes are plotted as a function of k,z in figure [5] (a). The collapse of all the mode shapes 
on the same function proves their self-similarity. As H. increases with Re, every experiment has a contribution to the 
lower region of the graph. There are consequently many more points at the bottom of the curve, which explains the 
higher scatter of the experimental data on this part of the plot. 

The evolution of the self-similarity coefficient k, with (. is plotted on figure [5] (c). It strongly depends on the string's 
reduced length (., and exhibits many discontinuities. Each one of the continuous parts of this indented curve was 
observed to correspond to a different mode of vibration of the hanging string. For instance, the two cases shown in 
figure [3] (a) and (&) correspond to the points A and B in figure [5] (c), where the second and fourth structure modes 
are excited, respectively. The discontinuities in the evolution of the self-similarity coefficient k therefore seem to be 
associated with changes in the modes of vibrations. Further insight on this mode selection is provided in the next 
section using a wake-oscillator model. 

V. A LINEAR WAKE-OSCILLATOR MODEL 

Wake-oscillator models are reduced-order models j24| that have been proven recently to render well the main features 
of the dynamics of a continuously flexible structure in VIV, in particular its frequency and mode shape [l0|. Moreover, 
the linear stability analysis of such a model has already shown a good ability to predict the dominant zone of each 
mode as well as its particular mode shape and frequency for tensioned cables in VIV [23] . This approach is adapted 
here to the case of a hanging string, where tension varies. 

The dimensionless equation for the cross-flow displacement of the hanging string reads as 

y-{zy')' ^ ffiu.d, (7) 

where f fluid stands for the fluid forces on the solid apart from the added mass effect, which is already included in the 
deflnition of the effective mass per unit length rut [23| • The fluid forcing is therefore proportional to a fluctuating lift 
coefficient, f fluid = Mq where q = 2Cl/Clo and M = Cio/lS/^Tr^S't^, Clo being the mean value of the fiuctuating 
lift coefficient acting on the string and /x = mt/pD'^ its mass ratio. The evolution of q satisfies a harmonic oscillator 
equation forced by the string's dynamics, 

q + q^ Ay. (8) 

Following Refs. [2l|, [2^, the values of the model parameters are taken as A = 12, Clo — 0.3. Looking for solutions 
of the form [y,q\ = Re[e'"* (iy9 (z) , f/; (z))] , it results from equation ([5]) that ip = AuP (f / {iJ^ — 1^ . Equation ([7]) 
consequently becomes 

(V)' + fcV = 0, (9) 

where the wavenumber k and the frequency uj are linked through the dispersion relation 

w4 + w2(AAf-l-fc2)+fc2=G. (10) 
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FIG. 5: Self-similarity of the mode shapes, (a) Superposition of all the experimental mode shapes as a function of the rescaled 
spanwise coordinate. (6) Comparison between the experimental mode shapes (grey) and the prediction of the linear model in 
section FVl (black continuous line), (c) Evolution of the self-similarity coefficient k with the dimensionless length of the hanging 
string I, obtained from the experimental data (symbols) and using the linear model (continuous line) 

The displacement of the string's lower end must remain finite, leading to the solution ip (z) = Jq (^h^fz) for equa- 
tion ([9|) , where Jq is the zero-order Bessel function of the first kind [1^, [23 . The finite length of the string and the 
boundary condition at the fixed upper end restricts k to discrete values, namely 
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an 
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where a„ are the successive zeros of the Bessel function Jq. 

For given values of n and i, the dispersion relation, equation (jlOp . provides the growth rate of an unstable mode as 



-y/-kl + 2kn + AM-l. 



(12) 



From equations (ITT]) and (fT^ , it can be found that mode n is the most unstable mode for values of the dimensioness 
length i in the interval 



1 2 1 2 

— (a„_i+a„) < £ < — (a„ + a„+i) . (13) 

The most unstable mode shape is Jq (2fc„y'z) with fc„ = a„/2v^ over the interval defined by equation (J13p . This 
form is consistent with the self-similarity of the mode shapes which has been observed experimentally in the previous 
section. The mode shapes are truncature of a unique function Jq {2kn\/z) and the self-similarity coefficient is linked 
with the wavenumbers A:„ by the following relation 



(14) 



^TCf 



where fc^e/ is the most unstable wavenumber corresponding to (/?re/- The value of fc^e/ = 1.017 is derived according to 
the preceeding analysis, equation (jlip. The most unstable linear mode is shown in figure [5] (6), in excellent agreement 
with the mode shapes obtained experimentally. Furthermore, the values of the self-similarity coefficient k obtained 
analytically, equation (ITU) , are also plotted in figure [S](c). Again, the agreement of the predictions obtained using the 
linear wake-oscillator model with experiments is excellent. The stability analysis explains and predicts analytically 
each feature of the dependence of k with €, namely the origin and location of the discontinuities (change in the most 
linearly unstable mode) and the variations of k. 

The real part of the most unstable frequency is also derived in the same manner as the growth rate but is not shown 
here for brevity. Its evolution with £, centered around a value close to 1, exhibits some discontinuities because of the 
mode selection, similarly to the wavenumber fc. The impact of the mode selection may explain part of the dispersion 
of the Strouhal number measurements in figure Ul 

VI. CONCLUSION 

VIV of flexible slender bodies are characterized by the possibility of successive excitation of higher structural modes 
as the flow velocity is increased. In this paper, the dynamics of a hanging string in VIV is investigated experimentally 
to explore the effect of non-uniform tension. Over the range of Rt considered, standing vortex-induced waves were 
observed, with a concentration of the highest amplitudes in the regions of lowest tension. Excitations of up to the 
7th mode of the string were observed for relatively low Re. The experiments prove the self-similarity of the spatial 
distribution of the vibrations and these dynamics were fully characterized. Furthermore, the variations of the self- 
similarity coefficient with the reduced length of the string was explained, using a linear stability analysis, as the result 
of the successive lock-ins of the different eigenmodes of the string. This standing wave structure breaks down at 
higher Re: travelling waves are indeed observed in the lower part of the chain that exhibits the largest amplitude of 
vibrations. This phenomenon, that is beyond the scope of the present paper, should be further studied to identify 
wether it is intimately linked to the system used in the present experiments or a fundamental property of hanging 
strings in high Re flows. 
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